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Abstract 

The issue of predicting solar flares is one of the most fundamental in physics, 
addressing issues of plasma physics, high-energy physics, and modelling of complex 
systems. It also poses societal consequences, with our ever-increasing need for accu- 
rate space weather forecasts. Solar flares arise naturally as a competition between 
an input (flux emergence and rearrangement) in the photosphere and an output 
(electrical current build up and resistive dissipation) in the corona. Although ini- 
tially localised, this redistribution affects neighbouring regions and an avalanche 
occurs resulting in large scale eruptions of plasma, particles, and magnetic field. 
As flares are powered from the stressed field rooted in the photosphere, a study 
of the photospheric magnetic complexity can be used to both predict activity and 
understand the physics of the magnetic field. The magnetic energy spectrum and 
multifractal spectrum are highlighted as two possible approaches to this. 
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1 Introduction 



Solar flares are among the most energetic events in the solar system and 
have intrigued generations of physicists. These events influence a panorama of 
physical systems, from the photosphere of the Sun, through the heliosphere 
and into geospace. Flares occur in active regions in the solar corona, volumes in 
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the solar atmospheric plasma characterised by increased emission at 1,000,000 
K or more and containing kilogauss magnetic fields. Active regions are be- 
lieved to be formed through the convective action of subsurface fluid motions 
pushing magnetic flux tubes through the photosphere. These active region 
flux tubes are jostled around by turbulent photospheric and sub-photospheric 
motions and, when conditions are right, the active region produces a flare. 
The energy of these events (typically ^ 10^^ erg) is thought to come from 
the energy stored in the magnetic field being suddenly converted into other 
forms, accelerating particles to near-relativistic speeds, and creating tempera- 
tures in excess of 10,000,000 K. The precise conditions required to create these 
enormously energetic events are as yet unknown. 

Solar flares are one of the main aspects of the larger fleld of space weather. 
Space weather generally refers to the interaction of solar particles and mag- 
netic flelds with the Earth's magnetosphere and upper atmosphere. Under- 
standing this interaction is of considerable practical importance because tech- 
nological systems, such as communications and navigation satellites, can suffer 
interruptions or permanent damage. Energy releases as extreme as solar flares 
and coronal mass ejections (CMEs), both in terms of their physics and po- 
tential human impact, are of fundamental importance in any consideration of 
space weather. This has long been recognized, from their first observation (R. 
C. Carrington, September 1st, 1859, recorded a white light fiare), to the incep- 
tion of sunspot classification systems (Hale 1919; Mcintosh 1990), to current 
space weather considerations (Baker et al. 2008). In this paper we concentrate 
solely on the solar flare aspect of space weather prediction. 

The concept of forecasting space weather is dependent on our ability to 
predict parameters of solar flares in advance using remote sensing data. The 
most important aspects of any solar fiare are size ( a proxy for energy re- 
leased in one flare), frequency ( a proxy for total energy released over some 
timespan), onset time, and location. An accurate onset time and location pre- 
diction would allow scientists to predict the geo-effectiveness of any event. 
Realistically a full accurate prediction is impossible with remote sensing data. 
However we can hope to obtain a probability measure of flare occurrence (at 
least we could predict an 'all-clear') by studying the magnetic fleld of active 
regions. Giovanelli (1939) showed that size, type (i.e., magnetic conflguration) 
and development of an active region could be used to predict flare activity. 
Kiinzel (1960) recognised the importance of the photospheric magnetic field 
measurements in the concept of adding the 6 configuration into the Mount 
Wilson classification. The presence of a 5 configuration, where large values of 
opposite polarity exist close together, was identified as a warning of the build 
up of magnetic energy stress, and hence the possibility of the occurrence of a 
large solar fiare. Mayfield & Lawrence (1985) showed a good correlation be- 
tween fiare energy and magnetic energy and found a 6 component doubled the 
probability of a large flare. Sammis, Tang & Zirin (2000) showed that large re- 
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gions classified as containing a S had a 40% chance of producing a large event 
of size XI or greater. These studies used this qualitative concept of magnetic 
complexity to show a maximum flare size could be estimated over some future 
period of time. An alternative 'Zurich' system (Keipenheuer, 1953), modified 
by Mcintosh (1990), has also shown some flare predictive quahties (Bornmann 
&; Shaw 1994) who show that even for the 60 different possible classes, flare 
predictions are of limited use for space weather purposes. 

There are two main approaches to the flare prediction problem. One is to 
combine existing active regions classifications systems and historical records. 
This approach is adopted by Sammis, Tang & Zirin (2000) to study the his- 
torical maximum flare size from the Mound Wilson classes. Gallagher et al. 
(2002) make a statistical estimate of flare occurrence above some threshold, 
based on a historical record of Mcintosh classes. Related to this, Bayesian 
studies (Wheatland 2001, 2004) of the waiting time distributions seem to do 
well at predicting the frequency of large events, although may over predict 
for more moderate size flares. Qahwaji & Colak (2007) and Colak & Qahwaji 
(2008) overcomes the qualitative nature of this approach somewhat by using 
machine learning to classify each active region and assign a flare prediction. 
Many of these approaches currently produce flare predictions in near-realtime 
which are made available online (www.swpc.noaa.gov, www.solarmonitor.org, 
spaceweather.inf.brad.ac.uk) . 

The second technique is to try to quantitatively assign an indicator of 
magnetic energy or complexity to active regions and test the evolution of this 
against the flares from the region. Abramenko et al. (2002, 2003) and Abra- 
menko (2005b) suggest the scaling index as a useful single parameter measure 
of magnetic complexity. McAteer, Gallagher & Ireland (2005a) suggest the 
fractal dimension seems to present a minimum threshold for flare size. Schri- 
jver (2007) suggest that an 'R- value', a weighted measure of flux near a neutral 
line, scales with solar flare size. Falconer, Moore, Gary (2002, 2003, 2006) pro- 
pose a number of different measures which scale with increased helicity and 
non-potentiality, and show how these can be used to predict CME produc- 
tivity. The role of helicity is also studied by Nindos & Andrews (2004) who 
connect the amount of helicity to whether a flare will be eruptive or confined. 
Leka & Barnes (2003) adapt a large number of measures of magnetic complex- 
ity and show, in a statistical sense, which measures may act as best predictors 
of flare occurrence. These studies are quantitative and mostly automated, so 
could conceivably provide near-realtime predictions. However the current lack 
of full-disc vector magnetic field information limits the usefulness of many of 
them for space weather prediction. 

In this paper we present the basic physics confirming the usefulness of 
using magnetic held information in solar flare prediction. We concentrate on 
characterizing active region magnetic field complexity in an attempt to begin 
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to understand which active region properties are important indicators of their 
activity. In Section [2] we discuss advantages and disadvantages of different 
types of magnetogram data. Section [3] contains a discussion of the conditions in 
the photospheric magnetic field and how this leads naturally to the algorithms 
presented in Section |H The algorithms of turbulence and fractals described in 
this work are based on original work by Kolmogorov (1941) and Mandelbrot 
(1983) and have since been found to ubiquitous in many areas of human and 
natural sciences (e.g., heartbeat dynamics, hydrology). These represent novel 
approaches to analyzing longitudinal magnetogram data and aim to generate 
physically motivated complexity measures. We conclude by discussing some of 
the drawbacks of these methods and possibilities for the future in Section [51 



2 Data 



In an ideal world, we could measure the full in-situ magnetic field vector 
from the photosphere to the corona, from which we could directly calculate 
any parameter of magnetic energy.. However, even in this seemingly perfect 
experiment, we would undoubtedly not find a perfect prediction of solar flare 
activity - this arises dues to the inherent non-linear nature of the solution 
of the Navier-Stokes equation, and dependence on small errors in the initial 
conditions. In our real world, we are limited to remote sensing of the magnetic 
field at one or two heights in the atmosphere, which brings about a number 
of extra issues of the resolution and coverage in the spectral, temporal and 
spatial domains (McAteer et al. 2005b). Spectrally, we would like to mea- 
sure the full vector field at multiple heights throughout the solar atmosphere. 
Temporally, we would like a full solar cycle of data from one instrument at a 
cadence higher than the evolution timescale of energy build-up. Spatially we 
would like to resolve the smallest features (i.e., elementary flux ropes), while 
covering the entire solar disk. Anything less than this ideal dataset and we 
have to be very careful of selection effects. Current instrumentation allows us 
to choose between full-disc longitudinal-measurements over a long timerange 
or partial disc vector field measurements of a few dozen regions. Specifically, 
the Michelson Doppler Imager (MDI; Scherrer et al. 1995) provides 96-min 
cadence, 1.96" /pixel, full disc images of the photospheric longitudinal (from 
Stokes V and I) field and has operated since 1997. Complementary to this, the 
Global Oscillations Network Group (GONG; Harvey et al. 1996) also provides 
full disc photospheric longitudinal magnetograms since 1995 and the Synop- 
tic Optical Long-term Investigations of the Sun (SOLIS; Keller et al. 2001) 
provides similar data for the chromosphere. Small field-of-view vector field 
measurements are available from the Imaging Vector Magnetograph (IVM; 
Mickey et al. 1996) and Hinode Solar Optical Telescope (SOT; Tsuneta et al 
2008). Soon we expect full disc vector magnetic field images from SOLIS, the 
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COronal Solar Magnetism Observatory (COSMO) and the Helioseismic and 
Magnetic Imager (HMI). 

To add to this quagmire of confusion, the choices of data and algorithms 
are intrinsically interlinked. For space weather purposes the data must be 
reliable, readily available, and easy to obtain. Several measurements of helicity 
and non-potentiality are tied to vector data and so we must wait on HMI 
before we can fully exploit their potential and then we must wait several years 
before having enough data to test any flare prediction system. In the meantime 
our algorithm of choice must be appropriate for the longitudinal data (e.g., 
MDI) while producing a physically meaningful measure of flare production. 
Furthermore, it must run efficiently and quickly. As the algorithms we discuss 
in Section m are naturally scale free, as the data are available in near- realtime, 
and as we find that a full-solar cycle of regular data is an essential aspect, we 
suggest MDI will provide us with the best data for fiare prediction over the 
next few years. 



3 The Magnetic Reynolds Number 



It seems that the magnetic field in active regions is the only viable means 
of storing and releasing the energy to drive solar fiares. A simple discussion of 
the equations used to describe this magnetic field leads to a couple of possible 
choices of characterising the field. 

The induction equation is given by eliminating the electric field between 
Faraday's law and Ohm's law 

^ = Vx(vxB)-V^, (1) 
ot a 

which further reduced by Ampere's Law to 

= V X (v X B) - V X (r/V X B), (2) 



where rj = l/fio<^ is the magnetic diffusivity. By Gauss' law, and by means 
of a simple vector identity (and assuming constant t]), this further reduces to 
the more familiar form of the induction equation, 

= V X (v X B) + r^V^B. (3) 

at 



This states that any local change in the magnetic field is due to a combination 
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of a convection and a diffusive term. The ratio of tfiese two terms is described 
by the magnetic Reynolds number, 

_ Vx{vxB) ^vl 



which acts as a indication of the couphng between the plasma flow and the 
magnetic field. For typical photosphcric values (/ ~ lO^Mm, v ~ 10ms^^,ri ^ 
lO^m^s"^), the magnetic Reynolds number is much greater than one. In this 
Rm ^ 1 regime, flux lines of the magnetic field are advected with the plasma 
flow, until such time that gradients are concentrated into short enough length 
scale that diffusion can balance convection (/ ^ fOOm), Ohmic dissipation 
becomes important, and magnetic reconnection can occur. Essentially the large 
Rm allows for the build up of energy, followed by a sudden release. Physical 
system with large Rm naturally lead to studies of turbulence and complexity 
as key parameters to developing a deeper understanding of the physics behind 
active region evolution and flare production. 



4 Turbulence and Complexity Measures 



The high magnetic Reynolds number in the photosphere suggests two pos- 
sible, interlinked, methods of quantifying complexity. These are studies of the 
scaling index (arising from fully-developed turbulence) and fractality (arising 
from self-similarity). Previous studies have shown that the local properties of 
the active region field - critical in many theories of activity - are lost in the 
common global definition of their diagnostics, in effect smoothing out varia- 
tions that occur on small spatial scales. While traditional Fourier and fractal 
methods have been used extensively in image processing, and are well un- 
derstood, they do not provide a complete diagnostic of the range of spatial 
frequencies or scales within an image. Hence measures that are sensitive to 
the small-scale nature of energy storage and release in the solar atmosphere 
are required. As fiows on the Sun exist in a state of fully developed turbu- 
lence, multiscale methods may be key to measuring and understanding the 
complex structures observed on the solar surface. In a similar fashion, multi- 
fractal, rather than monofractal, studies may be necessary to study the full 
complexity in locally inhomogeneous data. Importantly, these parameters can 
be related directly to predictions from theories of turbulence (Lawrence et al. 
1993) 

In this work we describe two complimentary techniques of multiscalar and 
multifractal methods. Both methods are based on the premise that physical 
systems cannot be adequately described by simple parameters, but require 
methods that capture their true complexity as a function of size scale. 
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4-1 Multifractal measures 



The fractal dimension of any object can be thought of as the self-similarity 
of an image across all scale sizes, or the scaling index of any length to area 
measure, 

A oc r, (5) 



where a is the singularity strength. However, a multifractal system will contain 
a spectrum of singularity strengths of different powers, 

A<xl^^''\ (6) 



and takes account of the measure at each point in space. Any measure distri- 
bution (e.g. magnetic field in an image) can be characterized by 

i^{q,T)=EJ2P^'e-^ (7) 

N 



where q, r can be any real numbers, and E is the expectation of the object 
consisting of N parts. In this form, ip is the coupled r -moment of the size e, and 
g-moment of the measure P. The three main multifractal indices commonly 
used to represent a non-uniform measure are then the : 

(i) Generalised Correlation dimensions (Grassberger & Procaccia 1983), Dq — 

(ii) singularity strength, a = dr/dq; 

(iii) Lcgcndre transformed f{ct) = qa — r. 

When applied to a traditional box-counting approach, it is useful to define the 
partition function, 

Z,{e) = Y.Piie), (8) 



such that T{q) — lime^Qlog{Z)/log{e), and any of the three representations 
above can be calculated. The terminology used in this description is deliber- 
ately similar to that use in turbulence studies (McAtccr ct al. 2007). The q 
moment plays the role of increasing the relative importance of the more in- 
tense parts of the measure as q is increased. In this way it acts as a microscope 
to investigate the different contributions made to the image at higher values 
of the measure. 

Another method of calculating the multifractal spectrum is based on the 
structure function (Parisi & Prisch 1985; Abramenko et al. 2002). This con- 
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sists of calculating the statistical moments of the field increments Sq{ T ^ J as a 
function of separation, r in order to determine the scaling exponents, 

Sq{r) ~ r^' (9) 



which are also directly related to the f{a) spectrum as (Muzy, Bacry & Ar- 
neodo 1993) 

D{h) = qh-C, + l (10) 



This formulation also provides a direct link to the Fourier scaling index (Sec- 
tion 113) as /? = 1 - Ce 

Another recent advance in this field uses the properties of the wavelet trans- 
form to order to overcome the limitations of the box-counting method (Conlon 
et al. 2009, Kestener et al. 2009). In this terminology two properties, D(/ig) and 
hq, are calculated separately by chaining the modulus maxima of the wavelet 
transform of an image, resulting in a D{h) spectrum which is directly related 
to the /(a) spectrum as 
(i) f{a) = Dihq), 
(n) hg = a- Edim, 

where Edim is the euclidian dimension {Edim = 2 for an image). This wavelet 
transform modulus maxima (WTMM) formulation offers many computational 
advantages over the box-counting and structure function approaches, includ- 
ing more accurate results at negative q and large positive q. 

Figure 1 shows a comparison of the Dg and /(a) spectra for a monofractal 
(top left; the Sierpinski carpet), a multifractal (middle), and a solar active re- 
gion magnetogram (right). The monofractal is described as a fiat Dg spectrum 
{Dg ~ 1.89 for all q)- it demonstrates the same complexity at all moments - 
and consequently a narrow /(a) spectrum. In the limit of an ideal algorithm, 
/(a) would consist of a single point at /(a = 2) = 2. The multifractal exhibits 
a monotonically decreasing Dg spectrum, which can be interpreted as meaning 
the image consists of series of fractals, each of which dominate at different mo- 
ments. The /(a) spectrum is wide and drops off rapidly which shows that the 
image contains a large number of singularity strengths each with a different 
power. Multifractals always show a decreasing Dg and more complex images 
contains more power at large singularity strengths. It is important to note 
that the degree of mutlifractality (e.g., the width and height of /(a)) reflects 
the range of fractals which exist in the image. 

The three segments of the MDI image all display a distinct drop off with 
increasing q, and hence are all multifractal. The complex segment shows a 
high Dg ~ 1.85 at large q, and a narrow, well-peaked f{a) spectrum. This in 
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Fig. 1. Top: An example of a monofractal (Sierpinski carpet; left), a multifractal 
(middle), and an aetive region (NOAA 9077) from MDI (right). Middle: The Dg 
and f(a) spectrum of the monofractal (crosses) and multifractal (diamonds). The 
solid lines correspond to the theoretical values. Bottom: The Dg and f{a) spectrum 
of three segments of the image of the active region 

is interpreted as meaning that the regions of strong field within this segment 
are highly complex, and contribute as much complexity to the image as the 
lower magnetic field parts. This can be contrasted with the plage segment of 
the data, which exhibits a Dg spectrum with a similar drop off at small q 
but continues to drop off at larger q. Correspondingly, the f{a) spectrum is 
much wider. Hence the strong field parts of the plage segment are rare and 
less complex. Finally, the quiet-Sun segment displays a Dg spectrum which is 
much lower at large q. The /(«) is very sparse at the high a values and drops 
off to very low /(a) at the lower a values. Hence the quiet Sun segment is 
mostly small, weak field. 

We can compare the differences in these spectra to the probabilities of 
flares occurring in each part of the image. Multifractal measures can be di- 
rectly compared to self-similar cascade or self-organized criticality models, 
both of which have been used to model energy release in solar flares (Lawrence 
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et al. 1993; Gcorgoulis ct al 2002). Fractal and multifractal algorithms have 
been apphed extensively to photospheric magnetic field data. Abramenko et 
al. (2002) found that the relative fraction of small scale fluctuation in the 
magnetic field contribute significantly more prior to fiaring. Abramenko et al. 
(2005a) found that active regions reach a critical state of intermittency prior 
to fiaring. McAtccr ct al. (2005a) tested a simple monofractal approach on the 
largest datasct to date, and found a minimum Dq = 1.2 is a necessary, but not 
sufficient requirement for M and X-class fiares. Conlon et al. (2008, 2009) find 
two distinct thresholds in the multifractal spectrum which correspond to the 
onset of fiare production in an otherwise quiet active region. They confirm a 
minimum Haussdorf dimension of D{h) = 1.2 is a necessary, but not sufficient 
requirement for M and X-class flares. Furthermore this must be accompanied 
by a minimum Holder exponent of -0.7. Physically this corresponds to a global 
restructuring of of the field distribution from Dirac like noise {h — —1) to a, 
step-function {h — 0) form. As the Holder exponent increases towards zero, 
salt-and-pcppcr type noise is replaced by the formation of gradients, permit- 
ting the build up of energy in the system 

4-2 Multiscalar measures 

It is recognised that small regions of fiux emergence / cancellation are vital 
in detailing the evolution of solar active regions. For this reason the wavelet 
transform can be adopted as it localised in space and hence allows for the 
detection of local image features. 

The continuous wavelet transform of an image, /(r) can be defined as 



where ^' is the mother wavelet, s is a term describing scale at a position 
r, and w{s,x) are the wavelet coefficients of the image. The mother wavelet 
can take several forms, depending on the application. Wavelet analysis retains 
the localized spatial information, providing vital information on the turbulent 
flow. As such, a wavelet analysis provides an indispensable complement to 
a multifractal analysis. The choice of mother wavelet is determined by the 
science requirements. A 'derivative of Gaussian' (also know as mexican hat) 
wavelet is useful for detecting sudden local changes in an image hence this can 
be use to detect the neutral line position and magnitude. A 'Haar' wavelet 
is useful in detecting and removing the intermittent component of an image 
(McAteer & Bloomfleld 2009), a vital tool in describing bursty behaviour. A 
discrete (e.g., 'a trous') wavelet can be used to efficiently detect fiux emergence 
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Fig. 2. A wavelet decomposition of the complex segment (region 1; top), plage seg- 
ment (region 2: middle), and quiet Sun segment (region 3: bottom) of the active 
region from Figure 1 at scales of 3.2 Mm (1st column), 9.6 Mm (2nd column), 
and 20.0 Mm (third column). The fourth column shows the energy spectrum as 
calculated from the Fourier technique (crosses) and the wavelet analysis (dots) 



as a function of size scale. 



Figure 2 shows the wavelet analysis decomposition (using the mexican hat 
wavelet) at three spatial scales for the same three segments of the MDI data 
as in Figure 1. The nine decompositions are all plotted in the same grayscale 
range. This makes it clear, that the energy (which is just the sum of the square 
of each value in the image) drops off from large k to small k for any image 
(reading across the rows). At any one scale, the energy is higher in the complex 
segment than in the plage, and is higher in the plage than in the quiet Sun 
(reading down the columns). The final column shows the energy spectrum as 
calculated using both the normal Fourier technique (adopted from Abramenko 
et al. (2005b)), and the wavelet technique (adapted from Hewett et al, (2008)). 
The calculation of the energy spectrum helps to characterize the intermittency 
of an image. Kolmogorov (1941) showed that the energy spectrum {E{k) of a 
system, scales as the wavenumber, k, 

E{k) - k-'^ (12) 
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where /3 = 5/3 for fully developed turbTilcnce. In Figure 2, /3 is calculated from 
a linear regression of the E{k) plot over the 3-lOMm range (as suggested by 
Abramenko et al. (2005b)). It is clear that the complex segment contains more 
power than the plage (and quiet Sun ) at all k and preferentially contains much 
more power at small k (hence large spatial scales) which results in a larger 
scaling index. Hence a large 3 is suggestive of increased complexity and more, 
larger, solar flares as demonstrated in Abramenko et al. (2005b). 

There are two main advantages of the wavelet technique over the Fourier 
technique, both of which are clear from Figure 2. Firstly, the linear range of 
the energy spectrum is much larger and hence the value of j3 calculated is 
much less dependent on identifying the correct range of linearity. The Fourier 
technique must exhibit a drop off for E{k) for small k (where the spatial scale 
is too large compared to the image size) and large k (where the spatial scale 
is not periodic across the image) . The inherent ability of the wavelet analysis 
to adapt to both small and large k makes the resulting energy spectrum more 
linear. Hewett et al. (2008) showed that the energy spectrum of an image can 
be accurately calculated over a larger range of scales using a wavelet transform. 
This work found a sudden onset of flares from an otherwise flare-quiet active 
region when the (3 suddenly dips from a previous value of —1 to a much steeper 
value of —3. The flare onset seems to occur when the spectrum passes thru 
the Kolmogorov index of —5/3. Furthermore the linear range of the scaling 
increases a few days prior to a flare. The linear range is theoretically known 
to increase with increasing Rm but the exact relationship is still under study. 

The second advantage is the retention of the spatial information in the 
the decomposition from which we can extract physical parameters at different 
scales. At the smallest spatial scale, there is a clear neutral hue near the center 
of the complex image. However at the largest scale, the neutral line towards 
the lower right of the image is more signiflcant. Ireland et al. (2008) use the 
wavelet transform to decompose magnetic field images into different scale sizes 
and show these are related to different Mount Wilson classes. They also find 
a significant difference in the distribution of gradients between flaring and 
non flaring regions and flnd this is maintained over all scales. The wavelet 
transform helps to encode the ideas of magnetic gradients and flux emergence 
over all length scales, not just those where the strongest fleld are found. 



5 Discussion 

We have discussed the necessity for accurate solar flare prediction and have 
reviewed the significant progress made over the past few years using param- 
eters of turbulence and complexity. It is clear that the more complex regions 
produce more larger flares, and by using multiscale and multifractal tech- 
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niques we may be able to capture this notion of complexity into one or two 
parameters. These take advantage of the excellent availability of a wide range 
of solar data both present and expected in the future. Data from new instru- 
ments will undoubtedly advance our understanding of active regions and the 
mechanisms through which they may affect human society. These new data 
will require new algorithms, e.g., calculating the multifractal spectrum of a 
vector field, which are currently being developed. This is expected to provide 
significant progress in predicting the true flare potential of individual active 
regions in near-realtime. 

There are a number of obvious problems with these techniques / data which 
must be addressed over the next few years. Firstly, we are only measuring 
the driving component of the build-up of energy - this does not inevitably 
lead to solar flares. The mechanism to initiate a flare can probably only be 
identified by studying the transverse component of the magnetic field in the 
corona. Secondly we are not fully aware of the natural timescale of the Sun. 
Although humans may prefer a 24-hour prediction, the Sun has a multitude 
of differing timescales interacting. It appears the driver of energy build-up 
may occur over days, with sudden inputs over hours being most significant. 
There may be a strong solar cycle dependence - predictive tools applicable to 
cycle 23 may not work for cycle 24. Furthermore fiares may not be Poisson 
distributed - it is well known that one of the best predictive signs of a large 
fiare is a previous large fiare. Thirdly we are not yet fully aware of what we 
mean by a large flare. We are only beginning to study the energy distribution 
(between radiation, particles, and the coronal mass ejection) in a flare event 
(Emslie et al. 2004). For space weather purposes, the geo-effective quality is 
a much more important statistic and this does not necessarily scale linearly 
with the scientists' measure (e.g., GOES class). We also have to be careful 
not to only study large events. A full physical understanding requires our 
models to work over all flare sizes. Lastly all algorithms need to be fully tested 
using climatological skill scores in order to make the step from scientiflcally 
interesting to space-weather prediction (Bloomfleld et al. 2009) 

The author wishes to thank the organisers of the sessions at European 
Space Weather 2008 for the invitation to present this work. This paper was 
greatly enhanced by discussions at the NASA all-clear workshop with Barnes, 
Leka, and Gourgoulis. The content of this paper was significantly improved by 
the comments of two anonymous referees. The authors thank the SOHO/MDI 
consortia for making data available. This research was supported by a grant 
from the Ulysses Ireland-France Exchange Scheme operated by the Royal Irish 
Academy and the Ministre des Affaires Etrangres. Paul A. Conlon is an IRC- 
SET Government of Ireland Scholar. R.T.James McAteer is the recipient of a 
Marie Curie Intra-European Fellowship under FP6. 



13 



References 

Abramenko, V.I., Yurchyshyn, V.B., Wang, H., Spirock, T.J., & Goode, P.R. 
Scaling behavior of structure functions of the longitudinal magnetic field in 
active regions on the Sun Astrophys. J., 557, 487-495, 2002. 

Abramenko, V.I., Yurchyshyn, V.B., Wang, H., Spirock, T.J., & Goode, P.R. 
Signature of an avalanche in solar flares as measured by photospheric mag- 
netic fields Astrophys. J., 597, 1135-1144, 2003 

Abramenko, V.I.. Multifractal analysis of solar magnetograms Sol. Phys., 228, 
29-42, 2005a 

Abramenko, V.I.. Relationship between magnetic power spectrum and flare 
productivity in solar active regions Astrophys. J., 629, 1141-1149, 2003 

Bornmann, P.L. & Shaw, D., Flare rates and the mcintosh active- region clas- 
sifications Sol. Phys., 150, 127-146, 1994 

Bloomfield, D.S., McAteer, R.T.J. , & Gallagher, P.T. A chmatological test of 
solar flare prediction algorithms Submitted to Astrophys., J. , 2009 

Baker, D.N., Balstad, R., Bodaue, J.M. et al. Severe space weather events- 
understanding societal and economic impacts: a workshop report National 
Research Council Report, 2008 

Colak, T.,. & Qahwaji, R., Automated Mclntosh-based classification of 
sunspot groups using MDI images Sol. Phys., 248, 277-296, 2008 

Conlon, P. A., Gallagher, P.T., McAteer, R.T.J. , ct al Multifractal properties 
of evolving active regions Sol. Phys., 248, 297-309, 2008 

Conlon, P.A., Kestener, P., McAteer, R.T.J., et al.. Characterising complexity 
in embedded systems, ii application to solar magnetic field measurements 
Submitted to Astronomy & Astrophysics, 2009 

Emslie, A.G., Kucharek, H., Dennis, B.R., et al Energy partition in two solar 
flarc/CME events Journal of Geophysical Research, 109, A18, 10104+, 2004 

Falconer, D.A., Moore, R.L, & Gary, G.A. Correlation of the coronal mass 
ejection productivity of solar active regions with measures of their global 
nonpotentiality from vector magnetograms: baseline results Astrophys. J., 
569, 1016-1025, 2002 

Falconer, D.A., Moore, R.L, & Gary, G.A. A measure from linc-of-sight mag- 
netograms for prediction of coronal mass ejections Journal of Geophysical 
Research, 108, 1380+, 2003 

Falconer, D.A., Moore, R.L, & Gary, G.A. Magnetic causes of solar coronal 
mass ejections: dominance of the free magnetic energy over the magnetic 
twist alone Astrophys. J., 644, 1258-1272, 2006 

GAL02] Gallagher, P.T., Moon, Y.-J.& Wang, H. Active-region monitoring 
and flare forecasting i. data processing and first results Sol. Phys., 209, 
171-183, 2002, 

Giovanelli, R.G., The relations between eruptions and sunspots. Astrophys. 

J., 89, 555-567, 1939 
Georgoulis, M.K., Rust, D.M., Bernasconi, P.N., & Schmieder, B., Statistics, 

morphology, and energetics of EUerman bombs Astrophys. J., 575, 506-528, 



14 



2002 

Grassbcrgcr, & Procaccia, Measuring the strangeness of strange attractors 
Physica D: Nonlinear Phenomena, 9, 189-208, 1983 

Hale, G.E., EUerman, P., Nicholson, S.B., & Joy, A.H. The magnetic polarity 
of sun-spots Astrophys J., 49, 153-178, 1919 

Harvey, J. W.; Hill, P.; Hubbard, R. P. The global oscillation network group 
(GONG) project Science, 272, 5266, 1284-1286, 1996 

Hewett, R.J., Gallagher, P. ., McAteer, R.T.J., et al Multiscale analysis of 
active region evolution Sol. Phys., 248, 311-322, 2008 

Ireland, J. Young, C.A., McAteer, R.T.J, et al Multiresolution analysis of 
active region magnetic structure and its correlation with the mount wilson 
classification and flaring activity Sol. Phys., 252, 121-137, 2008 

Keller, C.U. the SOLIS team The SOLIS Vector-Spectromagnetograph Ad- 
vanced solar polarimetry - theory, observation, and instrumentation, ASP 
Conference Proceedings, Vol. 236. Edited by Michael Sigwarth. San Fran- 
cisco: Astronomical Society of the Pacific 

Kestener, P., Conlon, P. A., McAteer, R.T.J. , et al Characterising complex- 
ity in embedded systems, i method and image segmentation Submitted to 
Astronomy & Astrophysics, 2009 

Kunzel, H. Die flare-haufigkeit in fleckengruppen unterschiedlicher klasse und 
magnetischer struktur Astron. Nachr., 285, 271-273, 1960 

Kiepenheuer, K.O. Solar activity The Sun, 322-466, Edited by Gerard P. 
Kuiper., Chicago, The University of Chicago Press, 1953 

Kolmogorov, A.N. Dissipation of energy in a locally isotropic turbulence Dok- 
lady Akad. Nauk SSSR, 32, 141, 1941 English translation in: American 
Mathematical Society Translations 1958, Series 2, Vol 8, p. 87, Providence 
R.l) 

Lawrence, J. K.; Ruzmaikin, A. A.; Cadavid, A. C. Multifractal measure of 
the solar magnetic field Astrophys J., 417, 805-811, 1993 

Leka, K.D., & Barnes, G. Photospheric magnetic field properties of flaring ver- 
sus flare-quiet active regions, i. data, general approach, and sample results 
Astrophys J., 595, 1277-1295, 2003 

Mandelbrot, B.B. Fractal geometry of nature The Fractal Geometry of Nature, 
New York, Freeman, 1983 

Mayfleld, E.B., Lawrence, J.K., The correlation of solar flare production with 
magnetic energy in active regions Sol. Phys., 96, 293-305, 1985 

McAteer, R.T.J. , Gallagher, P.T., & Ireland, J. Statistics of active region 
complexity: a large-scale fractal dimension survey Astrophys., J., 631, 628- 
635, 2005 

McAteer, R.T.J. , GaUagher, P.T., Ireland, J, & Young, C.A Automated 
boundary-extraction and region-growing techniques applied to solar mag- 
netograms Sol. Phys., 228, 55-66, 2005 

McAteer, R.T.J., Young, C.A., Gallagher, P.T., & Ireland, J. The bursty na- 
ture of solar flare x-ray emission Astrophys., J., 662, 691-700, 2007 

McAteer, R.T.J., & Bloomfleld, D.S. The extraction of shocks using a local 



15 



intermittency measure Submitted to Astrophys., J., 2009 
Mcintosh, P.S. The classification of sunspot groups Sol. Phys., 125, 251-267, 
1990 

Mickey, D. L., Canfield, R., LaBonte, B. J., Leka, K. D., Waterson, M. F., 
& Weber The Imaging Vector Magnetograph at Haleakala Sol. Phys., 168, 
229-250, 1996 

Muzy, J.F., Bacry, E., & Arneodo, A. Multifractal formahsm for fractal sig- 
nals: the structure-function approach versus the wavelet-transform modulus- 
maxima method Phys. Rev. E, 47, 875-884, 1993 

Nindos, A. & Andrews, M.D., The association of big flares and coronal mass 
ejections: what is the role of magnetic helicity? Astrophys., J., 616, L175- 
L178, 2004 

Parisi, G., & Frisch, U. Turbulence and predictability of geophysical fluid dy- 
namics and climate dynamics. Proceedings of Enrico Fermi Varenna Physics 
School, LXXXVIII, edited by M. Ghil, , North Holland, Amsterdam, Nether- 
lands, 1985 

Qahwaji, R,& Colak, T. Automatic short-term solar flare prediction using 
machine learning and sunspot associations Sol. Phys., 241, 195-211, 2007 

Sammis, I., Tang, F., & Zirin, H. The dependence of large flare occurrence on 
the magnetic structure of sunspots Astrophys., J., 540, 583-587, 2000 

Scherrer, P. H.; Bogart, R. S.; Bush, R. I. The Solar Oscillations Investigation 
- Michelson Doppler Imager Sol. Phys., 162, 129-188, 1995 

Schrijver, C.J., A characteristic magnetic field pattern associated with all ma- 
jor solar fiares and its use in fiare forecasting Astrophys. J. 655, L117-L120, 
2007 

Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. The Solar Optical Telescope 
for the Hinodc mission: an overview Sol. Phys. 249, 167-196, 2008 

Wheatland, S. Rates of fiaring in individual active regions Sol. Phys., 203, 
87-106, 2001 

Wheatland, S. A bayesian approach to solar flare prediction Astrophys., J., 
609, 1134-1139, 2004 



16 



